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We review the effects of clumping on the profiles of resonance doublets. By allowing the 
ratio of the doublet oscillator strcnghts to be a free parameter, we demonstrate that doublet 
profiles contain more information than is normally utilized. In clumped (or porous) winds, 
this ratio can lies between unity and the ratio of the /-values, and can change as a function 
of velocity and time, depending on the fraction of the stellar disk that is covered by material 
moving at a particular velocity at a given moment. Using these insights, we present the results 
of SEI modeling of a sample of B supergiants, C, Pup and a time series for a star whose terminal 
velocity is low enough to make the components of its Si VIAA1400 independent. These results 
are interpreted within the framework of the Oskinova et al. (2007) model, and demonstrate 
how the doublet profiles can be used to extract infromation about wind structure. 



1 Introduction 

Clumping has been apparent from time series for 
many years (Kaper et al. 1999, Prinja et al. 2002). 
In this contribution wc will first show how clumping 
affects doublet ratios and then use a time series to 
demonstrate the presence of clumping in winds. 



2 Clumping 

It is well known in the AGN community that dou- 
blet ratios are sensitive to clumping, with the ratio 
related to the covering factor of the clumped medium 
(e.g., Ganguly et al. 1999). For an extended source, 
the optical depth ratio determined from a doublet 
lies between the ratio of the blue and red /-values 
i a = /b//r) an d unity, depending on the covering 
factor of the source. Further, it is possible to in- 
terpret the observed doublet ratio within the frame- 
work of the Oskinova et al. (2007) model. In their 
model, the doublet components have an observed op- 
tical depth ratio of 
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which varies between 1 and a (see, Massa et al. 2003, 
for a definition of T ra d, and Oskinova et al. for def- 
initions of the other symbols). Consequently, if wc 
allow the ratio of /-values to be a free parameter, 



rad 



Ik 
fk 



R 
K eff 



(2) 



Thus, f' B ) ' f' R (determined by the fit) gives , 
and, hence tM . Oskinova et al. also relate the mea- 
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sured and smoothed opacities (or optical depths) 
and the clump optical depths as follows: K e ff/Kf — 



(1 



; )/tc- Since Mq (where M is the mass loss 



rate and q is the ionization fraction, see, e.g., Massa 
et al. 2003) should be derived from Kf, the observed 

Mq must be corrected as follows 



Mq = (Mq) obs 
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Fig. 1 shows how the "observed ratio" defines a point 
on the y-axis which gives Tq. then gives the ra- 
tio of the effective to actual opacity, nf^ 
to obtain K? f t <SC nf and M's near expected values, 



/k^. Note: 



requires t<j > 5, which implies k^^/k^j < 1.1. 




Figure 1: Ratio of effective opacities, 
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(top) and effective to smooth opacities, 

K eff/ K f 
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(bottom) versus clump optical 



depth, Tq , for a doublet with a = 2. 
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the actual ratio of /-values, 2.02. The solid fit var- 
ied the ratio, giving a best fit value of 1.84. This 10% 
change clearly improves the fit, making the blue ab- 
sorption weaker relative to the red. The fit gives 
Tq ~ 0.5 and nf^/n^ ~ 0.5. Thus, the correction 

to Mq is only 1.3, increasing ratio of P V to radio 
mass loss rates given by Fullerton et al. (2006) from 
0.11 to 0.14 - still far smaller than expected, even if 
q(P v) ~ 0.5. Thus the mass loss rate in this case is 
truly smaller than expected. 
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Figure 2: Derived ratios of oscillator strengths ver- 
sus T e ff for Si iv AA1400 in B supergiants 
with wind lines 0.3 < r ra d < 5 (weaker 
have inadequate optical depth information 
and stronger are too saturated). 
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Figure 3: P v in £ Pup. Points are data, dashed 
curve is SEI fit with fixed /-values, solid 
curve allowed the ratio of /-values to vary. 

It must now be determined whether this addi- 
tional information is present in the observed profiles. 
As a first test, we fit the Si IVAA1400 doublet in the 
sample of B supergiants given in Prinja et al. (2005) 
using a variable /-value. The results are shown in 
Figure 2. Notice that all of the /-value ratios lie 
between 1 and 2, as expected for a clumped wind. 

As a second test, we fit Copernicus data of the 
P VAA1117, 1128 doublet in £ Pup. Two least 
squares SEI (Lamers et al. 1987) fits are shown 
in Fig. 3. Both use a T eff = 40kK, logg = 3.5 
TLUSTY model photosphere. The dashed fit uses 
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Figure 4: Dynamic spectra of the observed (left) 
and modeled (right) spectra of the blue 
(top) and red (bottom) components of 
Si IVAA1400 in HD 47240. The spectra are 
normalized by their means. 
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3 Clumping and time series 

Time series for stars with Voo > 2cAX/X, where AA 
is the doublet separation, present a valuable test for 
the effects of clumping on wind lines. The Bl lb, 
HD 47240, is such a star with = 980 km s _1 . 
Consequently, each of its Si iv components are ef- 
fectively independent. Thus, they were fit indepen- 
dently. The result is equivalent to allowing fb/fr to 
be a function of velocity. The results are shown in 
Fig. 4, which shows the observed and modeled nor- 
malized dynamic spectra versus spectrum number 
for each component, and Fig. 5, which shows the un- 
normalizcd optical depths and their ratio. Several 
points are noteworthy. First, both components show 
similar, persistent velocity dependent structure. The 
fits of many B supcrgiants given by Prinja et al. 
(2005) had similar structure appearing in different 
ions in the same star. Thus, this structure may 
be real. Whether it results from density inhomo- 
geneties or velocity plateaus, cannot be determined. 
Second, when r ra< j is greater than a few tenths and 
well-defined, the ratio of optical depths varies be- 
tween 1 (clumped) and 0.5 (unclumped). Third, it 
appears that the ratio (clumping) decreases at ve- 
locities where the density decreases. Fourth, there 
seems to be a general decrease in the clumping at 
larger velocity, similar that seen in O stars by Puis 
ct al. (2006). 

4 Discussion 

It has been demonstrated that resonance doublets 
contain more information than is usually exploited, 



that this information is related to clumping, and 
that it can be interpreted by the Oskinova et al. 
(2007) model. Our ultimate goal is to apply these 
techniques to a large number of stars with a range 
of stellar and wind parameters. This will allow us 
to examine how these relate to the empirically de- 
termined clumping parameters and provide clues to 
the physical agents responsible for wind clumping. 
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